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In 1978, R. Apery [1], [2] has given sequences of rational approximations to ((2) and £(3) yielding 
the irrationality of each of these numbers. One of the key ingredient of Apery's proof are second-order 
difference equations with polynomial coefficients satisfied by numerators and denominators of the above 
approximations. Recently, V. N. Sorokin [3] and this author [4], [5] have got independently a similar second- 
order difference equation for C(4). This recursion does not give diophantine approximations to £(4) = 7r 4 /90 
proving its irrationality, however it presents an algorithm for fast calculation of this constant. The aim of 
this note is to highlight a possible generalization of the above results for the number C(5), the irrationality 
of which is not yet proved since nowadays. 

1. Statement of the main result. Consider the difference equation 

(n + lfa {n)q n+1 + oi(n)g„ - 4(2n - l)a 2 (n)g n _i - 4(n - l) 4 (2n - l)(2n - 3)a (n + l)q n - 2 = 0, (1) 
where 

a (n) = 41218n 3 - 48459n 2 + 20010™ - 2871, 

ai(n) = 2(48802112n 9 + 89030880n 8 + 36002654n 7 - 24317344n 6 - 19538418n 5 
+ 1311365n 4 + 3790503n 3 + 460056n 2 - 271701n - 60291), 

a 2 (n) = 3874492n 8 - 2617900n 7 - 3144314n 6 + 2947148n 5 + 647130n 4 
- 1182926n 3 + 115771n 2 + 170716n - 44541, 

and define its three linearly independent solutions {q n }, {Pn}, and {p n } by the initial data 

87 1190161 

Qo = -1, qi = 42, q 2 = -17934, p = 0, p x = —, p 2 = — — , 

101 _ 344923 

p = 0, Pl = —, p 2 = — — 

2 16 

(here and below a symbol {x n } means the sequence {x n }^ =0 = {xo,xi,x 2 , •••})• 
Theorem 1. The sequences 

= <?nC(5) ~Pn, tn = q n C{?>) ~ Pn, Tl = 0, 1, 2, . . . , (2) 
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that also satisfy the difference equation (1), are of constant sign: 

£ n >0, I n <0, n = l,2,... 
and the following limit relations hold: 

log|C | 



(3) 



lim 



n— s-oo n 



lim 1<3g — = log \fj, 2 \ 



1.08607936. 



n— »oo Tl 



log\q n \ log\p n \ r log\p n \ 
lim = lim = lim = log |/2 3 1 , 



n— >oo n 



n— >oo n 



n— >oo fi 



where 



(4) 



(5) 



//i = -0.02001512..., /i 2 = 0.33753726 ... , ^ = -2368.31752213... 

are i/ie roots of the characteristic polynomial fi 3 + 2368/t/ 2 — 752^ — 16 of the recursion (1). 
Note that the sequences {q n }, {p n }, and {p n } are of alternating sign: 

(-ir"V>0, (-1)"-V>0, (-l)«- 1 p n> 0, n = l,2,.... 

Theorem 1 gives an algorithm for fast calculation of the number C(5)- Namely, the sequence of rational 
numbers p n /q n converges to ((b) with speed l^/^s] < 1.42521964 • 10~ 4 (see the table). 



n 


Pn/Qn 


IC(5) -Pn/q n \ 








1.036927755... 


1 


29 
28 


0.001213469... 


2 


24289 
23424 


0.000000182... 


3 


7682021239 
7408444032 


< 2.80- 10" 11 


4 


24943788950905 
24055474286592 


< 4.13 • 10" 15 


5 


81875586674776013003 
78959779279372800000 


< 6.02 • 10" 19 


6 


282653756112686336975107 
272587704119854963200000 


< 8.71 • 10" 23 


7 


215903781003833520407770175189 
208214873150908926517286400000 


< 1.26 • 10" 26 


10 




< 3.71 • 10" 38 


20 




< 1.32 • 10" 76 


50 




< 5.52 • 10" 192 



Further construction of the above solutions of the difference equation (1) leads to the inclusions 

4D 2 n q n eZ, 4D 7 nPn eZ, 4D*p n eZ, n = l,2,... 

(see relations (14) and (15)), where D n denotes the least common multiple of numbers 1,2, ... ,n, while 
straightforward calculations on the basis of the recursion (1) show that 



q n £ Z, 2D a nPn G Z, 2D 3 n p n € Z, n = 1, 2, . . . . 



(6) 



The inclusions (6) do not allow oneself to prove the irrationality of the number C(5), hence we do not 
accent our attention on arithmetic properties of the sequences {q n }, {p n }, {Pn} and stress only on an 
applied meaning of Theorem 1. 
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2. Auxiliary recursions. The very-well-poised hypergeometric series 
r n = n! 4 



OO 

n 



rn ~~ n - tS ^ J n?=o(*+;) 6 ' 

are Q-linear forms in 1, C(3), C(5): 

r„ = n„C(5) + u>„C(3) - v n , r n = u„((5) + w n ((3) - v n , n = 0,1,2,.... 
An easy verification shows that 

r = C(5), n = 9C(5) + 33C(3) - 49, r 2 = 469C(5) + ^C(3) - 

ro = C(3), ri = 2C(5) + 12C(3) - f , f 2 = 552C(5) + 1764C(3) - 

z lb 

Applying the algorithm of creative telescoping [6, Chapter 6] to the series (7) in a manner of the works [4], 

[5], [7], we arrive at the difference equations 

n(n + l) 5 b {n - l)u n+1 - 2nb 1 (n)u n - fe 2 (n)u„_i + 2(n - l) 5 (2n - l)b (n)u n - 2 = 0, (8) 

where 

b (n) = -a (-n), 6i(n) = o 2 (-n), 6 2 (n) = -ai(-n), 

and 

n 3 (n + l) 3 6 (n - l)u n+1 - 2nb\{n)u n - 6 2 (n)ii n _i + 2n(n - l) 4 (2n - 3)6 (n)u n _ 2 = 0, (9) 

where 

bo(n) = 41218n 7 + 35648n 6 - 932n 5 - 13190n 4 - 5128n 3 + 811n 2 + 957n + 174, 

6i(n) = 3874492n 12 - 14084302n n + 12425954n 10 + 8641603n 9 - 15230839n 8 - 1369195n 7 

+ 8618417n 6 - 623249n 5 - 2785973n 4 + 308165n 3 + 495325n 2 - 40670n - 37632, 
b 2 (n) = 2(48802112n 13 - 201803328n 12 + 267014032n n - 69927236n 10 

- 95912858n 9 + 37524471n 8 + 30257812n 7 - 9523224n 6 - 8524312n 5 

+ 2138687n 4 + 1507490n 3 - 398634n 2 - 111012n + 33408). 

Note that each of the recursions (8) and (9) has the same characteristic polynomial A 3 — 188A 2 — 2368A + 4; 
its roots Ai, A 2 , A3 ordered in increasing order of their moduli are related to the roots (5) as follows: 

Hi = AiA 2 , jU 2 = A1A3, ^3 = A 2 A3. 

Theorem 2. The sequences of the linear forms {r n } and their coefficients {u n }, {w n }, {v n } satisfy the 
difference equation (8). In addition, the inequalities 



r n > 0, u n > 0, w n > 0, v n > 0, n=l,2, (10) 



and the limit relations 



lim — - ■ — = logAi, lim — - ■ = lim — — — — = lim — - ■ — = log A3 (11) 

n— s-oo n n— >oo n n— >oo n n— >oo n 



hold. 
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Theorem 3. The sequences of the linear forms {f n } and their coefficients {u n }, {w n }, {v n } satisfy the 
difference equation (9). In addition, the inequalities 

f„<0, u n > 0, w n >0, v n >0, n = l,2,..., (12) 

and i/ie limit relations 



lim — = logAi, lim — - = lira — - = lira — - = log A3 (13) 

ra— »oo 71 ra— >oo 77, n— >oo 77 n— >oo 77 



hold. 



We raeatioa that Vasilyev's result [8] aad our theorem [4] (see also [9]) on the coincidence of very- 
well-poised hypergeometric series and a suitable generalization of Beukers' integral [10] lead one to the 
inclusions 

2u n eZ, 2D 2 n w n G Z, 2D 5 n v n E Z, 

n = 1,2, ... . 14 

2tj„ e Z, 2D 2 n w n € Z, 2D^ n G Z, 
Finally, the required sequences (2) are determined by the formulae 

C = *n^n - w n r n = {u n w n - u n w n )C(5) - {w n v n - w n v n ) , 

77 = 0,1,2,..., 

( n = u n r n - u n r n = {u n w n - u n w n )((3) - (u n v n - u n v n ), 

hence 

q n = u n w n - u n w n , p n = w n v n - w n v n , p n = u n v n - U n V n , H = 0, 1, 2, . . . . (15) 

For proving the inequalities (3) and the limit relations (4), it remains to use the estimates (10), (12), 
relations (11), (13), and Poincare's theorem. 

3. Conclusion remarks. Recently, another algorithm for fast calculation of ((b) and ((3) (also not 
producing good enough diophantine approximations to these constants), based on the infinite matrix 
product 



n 



, n it 
' ~2(2n+ 1) 2n(2n+ 1) n 1 

;/ 5 








C(5)\ 


.(: 





C,3,) 


Vo 








„=i 2(2n + 1) An 2 

\ 1 

has been suggested by R. W. Gosper [11]. However the letters [11] do not contain any description of analytic 
and arithmetic properties of the rational approximations so constructed. Mention also an algorithm due 
to E. A. Karatsuba [12] for calculation of the Riemann zeta function at positive integers. 

The above scheme of Section 2 also allows one to construct a third-order recursion for simultaneous 
rational approximations to ("(2) and ((3). Namely, consider the difference equation 

(n + l) 4 a (n)q n+1 - a\{n)q n + 4(2n - l)a 2 (n)q n _ 1 - 4(n - l) 2 (2n - l)(2n - 3)a (n + l)q n -2 = 0, (16) 

where 

a (n) = 946n 2 - 731n + 153, 

ai(n) = 2(104060n 6 + 127710n 5 + 12788n 4 - 34525n 3 - 8482n 2 + 3298n + 1071), 
a 2 (n) = 3784n 5 - 1032n 4 - 1925n 3 + 853n 2 + 328n - 184, 

and define its three linearly independent solutions {q' n }, {p' n }, and {p' n } by the initial data 
g/=l, q[ = U, q' 2 = 978, p' = 0, p\ = 17, p' 2 = p' = 0, pi = 23, p' 2 = 



A THIRD-ORDER APERY-LIKE RECURSION FOR C(5) 



5 



Theorem 4. The sequences {q' n }, {p' n }> an d {P' n } (of positive sign) as well as the sequences 

C = <?;C(3)-K, I' n = q' n m-P'n, n = 0,l,2,..., 
satisfy the limit relations 

lim = lim l ^K\ = lQg = _ L31018925 ; 

n — >oc ft n— »oo fl 

y loglgnl t jgg \Pn\ v IK I , , , 

lim = lim = lim = log ^3 , 

n — >oo fl n — >oc fl n — >oo ri 

where 

fi lt2 = 0.07260980 . . . ± f0.25981363 . . . , /i 3 = 219.85478039 . . . 
are the roots of the characteristic polynomial /i 3 — 220/i 2 + 32/i — 16 of the recursion (16). 
This time, auxiliary recursions are satisfied by the hypergeometric series 



, 2 ^ n? = i(fc-j) ., , 2 ^ n" =0 (fc-j) _ ni9 



^n-=o(fc+i) 3 ' n ■ fcn?=o(*+i) a 

that are Q-linear forms in 1, C(2), and C(3)- 

Similarly to the case of the difference equation (1), the explicit calculations lead to the inclusions 

q ' n ez, D 3 nP ' n eZ, D 2 n p' n eZ, n=l,2,..., 

that are quite better than one can expect. 
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